Bessel process, Schramm-Loewner 
evolution, and Dyson model 

— Complex Analysis applied to 
Stochastic Processes and Statistical Mechanics — * 

Makoto Katori ^ 
24 March 2011 

Abstract 

Bessel process is defined as the radial part of tlie Brownian motion (BM) in 
the D-dimensional space, and is considered as a one-parameter family of one- 
dimensional diffusion processes indexed by D, BES^^). First we give a brief review 
of BES(^\ in which D is extended to be a continuous positive parameter. It is 
well-known that Dc = 2 is the critical dimension such that, when D > (resp. 
D < -De)) the process is transient (resp. recurrent). Bessel flow is a notion such 
that we ree ard BES(^) with a fixed D as a one-parameter family of initial value 
X > 0. There is another critical dimension Dc = 3/2 and, in the intermediate values 
of D, Dc < D < Dc, behavior of Bessel flow is highly nontrivial. The dimension 
D = 3 is special, since in addition to the aspect that BES^^^ is a radial part of the 
three-dimensional BM, it has another aspect as a conditional BM to stay positive. 

Two topics in probability theory and statistical mechanics, the Schramm-Loewner 
evolution (SLE) and the Dyson model {i.e., Dyson's BM model with parameter 
/3 = 2), are discussed. The SLE*^^) is introduced as a 'complexification' of Bessel 
flow on the upper-half complex-plane, which is indexed by D > 1. It is explained 
that the existence of two critical dimensions Dc and Dc for BES*^^^ makes SLE^^^ 
have three phases; when D > Dc the SLE*^^^ path is simple, when Dc < D < Dc 
it is self-intersecting but not dense, and when 1 < D < Dc it is space-filling. The 
Dyson model is introduced as a multivariate extension of BES^^^. By 'inheritance' 
from BES^'^\ the Dyson model has two aspects; (i) as an eigenvalue process of a 
Hermitian- matrix- valued BM, and (ii) as a system of BMs conditioned never to col- 
lide with each other, which we simply call the noncolliding BM. The noncolliding BM 
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is constructed as a harmonic transform of absorbing BM in the Weyl chamber of type 
A, and as a complexification of this construction, the complex BM representation 
is proposed for the Dyson model. Determinantal expressions for spatio-temporal 
correlation functions with the asymmetric correlation kernel of Eynard-Mehta type 
are direct consequence of this representation. In summary, 'parenthood' of BES*^-^) 
and SLE(^\ and that of BES^^^ and the Dyson model are clarified. 

Other related topics concerning extreme value distributions of noncolliding dif- 
fusion processes, statistics of characteristic polynomials of random matrices, and 
scaling limit of Fomin's determinant for loop-erased random walks are also given. 

We note that the name of Bessel process is due to the special function called 
the modified Bessel function, SLE is a stochastic time-evolution of complex analytic 
function (conformal transformation), and the Weierstrass canonical product repre- 
sentation of entire functions plays an important role for the Dyson model. Complex 
analysis is effectively applied to study stochastic processes of interacting particles 
and statistical mechanics models exhibiting critical phenomena and fractal struc- 
tures in equilibrium and nonequilibrium states. 

Keywords Complexification, Multivariate extension, Conformal transformation. 
Random matrices. Entire functions 
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1 Family of Bessel processes 



1.1 One-dimensional and D- dimensional Brownian motions 

We consider motion of a Brownian particle in one dimensional space R starting from 
X e M at time t = 0. At each time t > particle position is randomly distributed, and 
each realization of path is labeled by a parameter u. Let Q be the sample path space 
and B^{t,uj) denote the position of the Brownian particle at time t > 0, whose path is 
realized a.s cu E fl. Let {fl, J-", P) be a probability space. The one- dimensional standard 
Brownian motion (BM), {B^{t,uj)}te[Q^oo),x G M, has the following three properties. 

1. B^{0,w) — X with probability one (abbr. w.p.l). 

2. For any fixed a; e B^{t,uj) is a real continuous function of t. In other words, 
B^{t) has a continuous path. 

3. For any sequence of times, to = < ti < ■ ■ ■ < tM,M e N = {1,2,3,...}, 
the increments {B^{ti) — B^{ti^i)}^^ are independent, and distribution of each 
increment is normal with mean m — and variance — ti — ti-i. It means that 
for any 1 < i < M and £ < r, 

P{B^{t,) - B^{t,_,) e [£,r]) = f pt^_t^_Mm, 

Ji 

where we define for a, 6 e IR 

Vt{h\a) = < V27rt (1-1) 
I 5{a-b), fori = 0. 



If we write the conditional probability as P(-|C), where C denotes the condition, the 
third property given above implies that for any < s <t 

F{B%t) e = a)= [ pt-s{h\a)dh 

J A 

holds C M,^ a e M. Then the probabihty that the BM is observed in a region C R 
at time ti for each i = 1, 2, . . . , M is given by 

P{B%ti) eAi,i^ 1,2,..., M)^ f dxi--- f dxMY[pu-u.M\xi-i), (1-2) 

JAi J Am j=i 

where xq = x. The formula (1.2) means that for any fixed s > 0, under the condition that 
B^{s) is given, {B^(t) : t < s} and {B^{t) : t > s} are independent. This independence 
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of the events in the future and those in the past is called Markov property ^. The integral 
kernel Pt{y\x) is called the transition probability density function of the BM. As defined 
by (1.1), Ptiy\x) is nothing but the probability density function of the normal distribution 
(the Gaussian distribution) with mean m = x and variance o"^ = t. It should be noted that 
Pt{y\x) = Pt{Av) "^^ti^) = PtiuW) is a unique solution of the heat equation {diffusion 
equation) 

d 1 

with the initial condition Uo{x) — S{x — y). Therefore, Pt{y\x) is also called the heat 
kernel. 

Let D E N denote the spatial dimension. For D > 2, the D-dimensional BM in 
M.^ starting from the position x — {xi, . . . ,xd) £ is defined by the following D- 
dimensional vector- valued BM, 

B^(t) = (i?r(t),i?r(t),...,i?^-(t)), te[0,oo), (1.4) 

where {B^'-{t)}^i are independent one-dimensional standard BMs. 



1.2 D-dimensional Bessel process 

For D E N, the D-dimensional Bessel process is defined as the absolute value (i.e., the 
radial coordinate) of the D-dimensional BM, 

X'^it) = \B^it)\ 

= ^sr(t)2 + ... + s^-(t)2, te[0,oc), (1.5) 

where the initial value is given by X^(0) = x = \x\ = \J x\ + ■ ■ ■ + > 0. By definition 
X^(t) is nonnegative, X^(t) e M+ = {x e M : a; > 0}. See Fig. 1. 

By this definition, X^ii) is a functional of D-tuples of stochastic processes {B^^ (^)}iii- 
In order to describe the statistics of a function of several random variables, we have to see 
'propagation of error'. For stochastic processes, by ltd's formula we can readily obtain 
an equation for the stochastic process that is defined as a functional of several stochastic 
processes. In the present case, we have the following equation, 

dX^(t)^dB{t) + ^^^y ie[0,oo), x>0. (1.6) 

The first term of the RHS, dB{t), denotes the infinitesimal increment of a one-dimensional 
standard BM starting from the origin at time t = 0, B{t) = B'^{t). It should be noted 

'^A positive random variable r is ealled Markov time if the event {r < u} is determined by the behavior 
of the process until time u and independent of that after u. The Brownian motion satisfies the property 
obtained by changing any deterministic time s > into any Markov time r in the definition of Markov 
property given here. It is called a strong Markov property. A stochastic process which is strong Markov 
and has a continuous path almost surely is called a diffusion process. 
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Figure 1: The D-dimensional Bessel process X^'(t) is defined as the radial part of the BM in 
the D-dimensional space. The initial value x of the Bessel process is the distance between the 
origin and the position from which the BM is started. 

that B{t) is a different BM from any B^'{t), 1 < i < D, which were used to define X^{t) 
in Eq.(1.5). Here we assume X^{t) > 0. Then, if D > 1, for an infinitesimal increment 
of time dt > 0, the second term in the RHS of (1.6) is positive. It means that there is a 
drift to increase the value of X^{t). This drift term is increasing in D and decreasing in 
X^{t). Since as X^{t) \ 0, the drift term oo, it seems that a 'repulsive force' is acting 
to the D-dimensional BM, B^(t), \x\ > to keep the distance from the origin positive, 
X^{t) = \B^(t)\ > 0, in order to avoid collision of the Brownian particle at the origin. 
Such a differential equation as (1.6), which involves random fluctuation term and drift 
term is called a stochastic differential equation (SDE). 

What is the origin of the repulsive force between the D-dimensional BM and the 
origin? Why B^(t) starting from a point x does not want to return to the origin ? 
Why the strength of the outward drift is increasing in the dimension D > 1 7 

There is no positive reason for B^{t) to avoid visiting the origin, since by defini- 
tion (1.4) all components B^^{t) enjoy independent BMs. As the dimension of space D 
increases, however, the possibility not to visit the origin (or the fixed special point) in- 
creases, since among D directions in the space only one direction is toward the origin (or 
the fixed special point) and other D — 1 directions are orthogonal to it. If one know the 
second law of thermodynamics, which is also called the law of increasing entropy, one will 
understand that we would like to say here that the repulsive force acting from the origin 
to the Bessel process is an 'entropy force'. (Note that the physical dimension of entropy 
[J/K] is different from that of force [J/m].) Anyway, the important fact is that, while the 
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variance (quadratic variation) of the standard BM is fixed as {dB{t)Y — dt for a given 
dt > 0, the strength of repulsive drift is increasing in D. Then, the return probabihty of 

X^{t), ,x > to the origin should be a decreasing function of D. 

Let Pt^\y\x) be the transition probability density of the D-dimensional Bessel process. 
We can show that, for any y e M+, u[^\x) = p[^\y\x) , x > solves the following partial 
differential equation (PDE) 

under the initial condition Uq^\x) = S(x — y), which is called the backward Kolmogorov 
equation for the D-dimensional Bessel process. We can see clear correspondence between 
the SDE (1.6) and the PDE (1.7). As shown by (1.3), the BM term, dB{t), in (1.6) is 
mapped to the diffusion term {\/2)d'^u\^\x) / dx^ in (1.7). In (1.7) the drift term is given 
by using the spatial derivative d/dx representing the outward drift with the coefficient 
{D — 1) /2x corresponding to the factor (D — l)/ {2X^{t)) of the second term in (1.6). The 
solution is given by 

' i^e-(^^+^^)/^%(^), t>0,x>0,y>0, 
P?\y\x) = I 2-t-+^r(z.+ l) ""' t>0,x^O,y>0 (1.8) 

. ^{y -x), t^O,x,y> 0, 

where Iiy{z) is the modified Bessel function of the first kind defined by 

W-E r(„ + i)r(„ + i + .) id (1-^' 



n=0 



with the gamma function r{z) — e ^u^ ^du, and the index u is specified by the 
dimension D as 

u^^—^ <^ D = 2(i/ + 1). (1.10) 

This fact that Vt^^ {y\x) is expressed by using lyi^z) gives the reason why the process X^{t) 
is called the Bessel process. 

When D = 3, i/=l/2by (1.10), and we can use the equality /i/2(-z) = v^7^ ^i^^ ^ — 
(e^ - e-^)/V27fz. Then (1.8) gives 

v't\y\x) = ^[pt{y\x) -Pt{y\ - a;)} (i.ii) 

for t > 0,x > 0,y > 0, where Pt{y\x) is the transition probability density (1.1) of BM. If 
we put pf'^iylx) = pt{y\x) — Pt{y\ — x), we see pf^^{0\x) = for any a: > 0, since BM is a 
symmetric process. 



6 



A 



A 



y 



-X 



X 




-y 



B 



(a) 



(b) 



Figure 2: (a) One realization of Brownian path from x > to y > is drawn (path A), which 
visits the nonpositive region M_ = {x G M : x < 0}. The path B is a mirror image of path A 
with respect to the origin x = 0, which is running from —x < to —y < 0. The first time when 
the path B hits the origin is denoted by r. The path C is a combination of a part of path B up 
to time r and a part of path A after r such that it runs from —x < to y > 0. There estabhshes 
bijection between path A and path C, which have the same weight as Brownian paths. Since 
the Brownian path A contributes to pt{y\x) and the Brownian path B does to pt{y\ — x), such 
a path from x > to y > visiting M_ is cancelled in pf^^{y\x). (b) Each path, which does 
not visit M_, gives positive contribution to pf^^{y\x), since in this case the weight of Brownian 
path A contributing to Pt{y\x) is bigger than that of Brownian path C contributing pt{y\ — x). 
In summary, pf^^{y\x) = Pt{y\x) —pt{y\ — x) gives the total weight of Brownian paths, which do 
not hit the origin and thus are not absorbed at the wall at the origin. 




Figure 3: The one-dimensional BM visits the origin frequently. For a Brownian path starting 
from X > (path A), its mirror image with respect to the origin is drawn (path B), which starts 
from —X < 0. If we observe the motion only in the nonnegative region M_|_, the superposition of 
Brownian paths A and B gives a path of a reflecting BM, where a reflecting wall is put at the 
origin. 
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As shown by Fig. 2, pf'^{y\x),x,y > 0, gives the transition probabihty density of the 
absorbing BM, in which an absorbing wall is put at the origin and, if the Brownian particle 

starting from a; > arrives at the origin, it is absorbed there and the motion is stopped. 
By absorption, the total mass of paths from a; > to y > is then reduced. The factor 
y/x appearing in (1.11) is for renormalization so that J^^pl'^\y\x)dy = l,Vt > 0,Va; > 0. 

We regard this renormalization procedure from pf^^ to p[^^ as a transformation. Since x 
is a one-dimensional harmonic function in a rather trivial sense A^^^x = (Px/dx^ = 0, we 
say that the three-dimensional Bessel process is an harmonic transform (/i-transform) of 
the one-dimensional absorbing BM in the sense of Doob [14]. This implies the equivalence 
between the three-dimensional Bessel process and 'the one-dimensional BM conditioned to 
stay positive\ We will discuss such equivalence of processes in Section 2 more detail. Here 
we put emphasize the fact that Pt^\o\x) = 0, Vx > 0. It means that the three-dimensional 
Bessel process does not visit the origin. When D — 3, the outward drift is strong enough 
to avoid any visit to the origin. Moreover, we can prove that for any x > 0, X^{t) — >■ oo 
as t — )■ cxD w.p.l and we say the process is transient. 

When D = 1, u = —1/2 by (1.10) and we use the equality I_i/2{z) = ^/qJttz cosh, z = 
(e^ -|- e~^)/\/27rz. In this case (1.8) gives 

p?{y\x) = pt{y\x) +pt{y\ - x) (1.12) 

for t > 0,x,y > 0. As shown by Fig.3, (1-12) means the equivalence between the one- 
dimensional Bessel process and 'the one-dimensional BM with a reflecting wall at the 
origin'. This is of course a direct consequence of the definition of Bessel process (1.5), 
since it gives X^{t) — \B^{t)\ in D = 1. The important fact is that the one-dimensional 
BM starting from x ^ visits the origin frequently and we say that the one- dimensional 
Bessel process is recurrent. (Remark that in Eqs. (1.6) and (1.7), the drift terms vanish 
when D — 1. So we have to assume the reflecting boundary condition at the origin when 
we discuss the one-dimensional Bessel process instead of the one-dimensional BM.) 

Now the following question is addressed: At which dimension the Bessel process 
changes its property from recurrent to transient ? 

Before answering this question, here we would like to extend the setting of the question. 
Originally, the Bessel process was deflned by (1.5) for D e N. We flnd that, however, 
the modified Bessel function (1.9) is an analytic function of u for all values of u. So we 
will be able to define the Bessel process for any positive value of dimension D > as the 
diffusion process in R_|_ such that the transition probability density function is given by 
(1.8), where the index u > —1 is determined by (1-10) for each value of D. (In the SDE, 
(1.6), we assume the refiecting boundary condition at the origin for < D < 2.) Now we 
introduce an abbreviation BES*^-^) for the D-dimensional Bessel process, D > 

^Another characterization of BES^^^ for fractional dimensions D is given by the following. Let e M 
and consider a BM with a constant drift v, By{t) + vt, which starts from y € R at time t = {). The 
geometric 5M with drift u is defined as exp(i?^(t) -\-vt),t > 0. For each t > 0, if we define random time 
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For BES^^^ starting from a; > 0, denote its first visiting time at the origin by 



= mi{t > : X'-'it) = 0}. 



(1.13) 



The answer of the above question is given by the following theorem. 



Theorem 1.1 (i 

(ii) D>2 = 



■) 



D>2 =^ T"" = oo^ x>Q, w.p.l. 
limX^(t) = oo/x > 0, w.p.l, i.e. the process is transient. 



(iii) 



D = 2 =^ infX^(t) = O/x > 0, w.p.iZ? 

t>o 

That is, BES^"^^ starting from x > Q does not visit the origin, but it can visit any 
neighbor of the origin. 

D <2 =^ < oo, '^x > 0, w.p.l, i.e. the process is recurrent. 



IV 



1.3 Bessel flow and Cardy's formula 

In the previous subsection, we have defined the BES*^^^ for positive continuous values of 
dimension D > and studied dependence of the probabihty law of process on D. Theorem 
1.1 states that the two-dimension is a critical dimension, 



for competition between the two effects acting the Bessel process, the 'random force' (the 
martingale term) and the 'entropy force' (the drift term) in (1.6) and (1-7): when D > Dc, 
the latter dominates the former and the process becomes transient, and when D < Dc, 
the former is relevant and recurrence to the origin of the process is realized frequently. 
Here we show that there is another critical dimension. 




In order to characterize the transition at we have to investigate dependence of the 

behavior of X^{t) on initial value x > 0. We call the one-parameter family {X^{t)}x>o 
the Bessel flow for each fixed D > 0. 

change t At hy 



where X^{At) is the BES*^^^ with D = 2{i' + 1) at time At starting from x = e^. The above formula is 
cahed Lamperti's relation [59, 88] 



D, = 2 




then the following relation is established, 



X''{At) = exp{B^ + 1^1), t>0, 
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time 



Figure 4: In the intermediate fractional dimensions, 3/2 < D < 2, there is a positive probability 
that two Bessel processes starting from different initial positions, < x < y < oo, return to the 
origin simultaneously, = T^. 

For < X < y, we trace the motions of two BES^'^^'s starting from x and y by solving 
(1.6) using the common BM, B{t),t > 0, 

D - 1 /•* ds 



D ~i r 

X^{t)=x + B{t) + ^- 

^'i^)-y^m + ^[^y t>o. 





By considering the couphng of the two processes, we can show that 
x<y =^ < Xy{t), t < w.p.l 

=^ T^' < ry w.p.l. 

The interesting fact is that in the intermediate fractional dimensions, < D < Dc, 
it is possible to see the coincide = even for x < y. See Fig.4. 

Theorem 1.2 For < x < y < oo, 

(i) D<3/2 =^ T"" <Ty w.p.l. 

(ii) 3/2 < < 2 =^ P(T^ = Ty) > 0. 



Theorem 1.2 (ii) is obtained by proving that, for < x < y < cxd,3/2 < D < 2, the 
event 

xy{t)-X''{t) 

sup ^y,.,. < oo (1.14) 
10 



occurs with positive probability. If (1.14) holds, < oo s.t. X^(t) < (1 + c)X^(t),0 < 
t < and thus X^(i) = ^ Xy{t) = ^ = T^^. We can confirm that the 
difference between the event {T^ — T^} and the event (1-14) has probability zero (see 
Section 1.10 in [60]). 

A striking fact is the following exact formula: for3/2<L'<2,0<a;<|/<oo, 



* ' mD-l))T(2-D)\ y ) 



2D-3 



xF{2D-3,D-l,2{D-iy,y—^], (1.15) 



where F{a, (3, 7; z) is Gauss' hypergeometric function 



i=0 



with the Pochhammer symbol {c)i = c(c + 1) • • • (c + i — 1). If we set D — 5/3, (1.15) 
becomes a version of exact formula for a physical quantity called 'crossing probability' that 
Cardy derived in a critical percolation model [10, 11]. Cardy's formula has been extended 
in the context of SLE (see Section 6.7 of [60]), but I think that this exact formula for the 
Bessel flow can be also called Cardy 's formula. 



1.4 Schramm-Loewner evolution (SLE) as complexification of 
Bessel flow 

Now we consider an extension of the Bessel flow X^{t) deflned on M_|_ to flow on the 
upper-half complex-plane M.— {z — x + \f—\y : x e M, |/ > 0} and its boundary dW — R. 
We set Z\t) = X\t) + ^f^Y^t) e H \ {0} = H U R \ {0}, t > and complexificate 
(1.6) as 

dZ^(t) = dS(i) + ie[0,oo) (1.16) 

with the initial condition 

Z^(0) = ^ = x + e H \ {0}. 

The crucial point of this complexiflcation of the Bessel flow is that the BM remains real. 
Then, there is asymmetry between the real part and the imaginary part of the flow in H, 



11 




Assume D > 1. Then as indicated by the minus sign in the RHS of (1.17), the flow is 
downward in H. If the flow goes down and arrives at the real axis, the imaginary part 
vanishes, Y^{t) = 0, and Eq.(1.17) is reduced to be the same equation as Eq.(1.6) for the 
BES(^), which is now considered for M \ {0}. If D > 2, by Theorem 1.1 (ii), the flow on 
M\{0} is asymptotically outward, — )■ ±oo as t — )■ oo. Therefore, the flow on H will be 
described as shown by Fig.5. The behavior of flow should be, however, more complicated 
when Dc = 3/2 < D < Dc and 1 < D < D^. 
For z G H \ {0},t > 0, let ^ 

gt{z) = Z'{t)-B{t). (1.19) 

Then, Eq.(1.16) is written as follows: 

d , , D-1 1 , , 

-g,{z) = — TY—^^ ^>0, 1.20 

at 2 gt{z) + B[t) 

with the initial condition go^z) = z eM. \ {0}. For each 2; G H \ {0}, set 

= inf{t > : gt{z) + B{t) = 0}, (1.21) 

then the solution of Eq.(1.20) exists up to time T^. For t > we put 

Ht = {zeM:T' > t}. (1.22) 



■^Since B'^{t),X^{t) and Z'^{t) are stochastic processes, they are considered as functions of time t>0, 
where the initial values x and z are put as superscripts {B{t) = B^{t)). On the other hand, as explained 
below, gt is considered as a conformal transformation from a domain TJt C H to H, and thus it is described 
as a function of z G Ht\ gtiz), where time i is a parameter and put as a subscript. 
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This ordinary differential equation (1.20) involving the BM is nothing but the cele- 
brated Schramm- Loewner evolution (SLE) [77, 60]. For each t > 0, the solution gt{z) of 
(1.20) gives a unique conformal transformation from Hf to H: 



gt{z) : Ht ^M, conformal 



such that 



gt{z) ^ z + 



a{t) 



z 




z ^ oo 



with 





z->0,zeH 



and using basic properties of BM, we can prove that 7 = 7[0, 00) = {7(t) : t G [0, 00)} G H 
is a continuous path w.p.l running from 7(0) = to 7(00) = 00 [74]. The path 7 obtained 
from the SLE with the parameter D > 1 is called the SLE^^^ path ^. 

The dependence on D of the Bessel flow given by Theorems 1.1 and 1.2 is mapped to 
the feature of the SLE^^^ paths such that there are three phases. 

[phase 1] When D > = 2, the SLE^^) path is a simple curve, i.e., 7(5) 7^ 7(t) for 
any < s 7^ t < 00, and 7(0, 00) G H (i.e., 7(0, 00) Pi M = 0). In this phase. 



For each t > 0, gt gives a map, which conformally erases a simple curve 7(0, t] from 
H, and the image of the tip j{t) of the SLE path is —B{t) G M = dM. as given 
by (1.23). As shown by Fig.6, it implies that the 'SLE flow' in H is downward in 
the vertical (imaginary- axis) direction and outward from the position —B{t) in the 
horizontal (real-axis) direction. Since Z''^{t) = gt{z) + B(t) by (1.19), if we shift 
this figure by B{t), we will have the similar picture to Fig.5 for the complexificated 
version of Bessel flow for D > 2. 

^Usual parameters used for the SLE are k = A/{D - 1) [77] or a = (£» - l)/2 = 2/k [60]. If we set 
9t{z) = \/Kgt{z) and B{t) = —W{t) in (1.20), we have the equation in the form [77], 



Note that ^JllW{t) has the same distribution with B{Kt), which is a time change t nt oi the one- 
dimensional standard BM. In SLE, the Loewner chain for H is driven by a one-dimensional BM, which is 
speeded up (or slowed down) by factor k compared with the standard one. (The parameter k is regarded 
as the diffusion constant.) 



Ht^m\-f{o,t], t>o. 




gt{z)-^W{ty 



2 
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7(0 




Figure 6: (a) When D > 2, the SLE^^) path is simple, (b) By gt, the SLE^^) path is erased from 
H. The tip of the SLE^^) path 7(t) is mapped to gtij{t)) = -B{t) G M. The flow associated by 
thus conformal transformation is shown by arrows. 

[phase 2] When = 3/2 < D < = 2, the SLE^^) path can osculate the real axis, 
P(7(0,t] n M 7^ 0) > 0, Vt > 0. Fig.7 (a) illustrates the moment t > such that the 
tip of SLE'^'^^ path just osculates the real axis. The closed region encircled by the 
path 7(0, t) and the line [7(t),0] G M is called an SLE hull at time t and denoted 
by Kt. In this phase 

Ht = m\Kt, t>0. 

That is, gt{z) is a map which erases conformally the SLE hull from H. We can think 
that by this transformation all the points in Kt are simultaneously mapped to a 
single point —B{t) G M, which is the image of the tip 7(t). (We say that the hull Kt 
is swallowed. See Fig.7 (b).) By definition (1.22), the moment when Kt = M.\Ht 
is swallowed is the time at which the equality Z^{t) = gt{z) + B{t) =0 holds 
'^z G Kf. (Then the RHS of (1.20) diverges and all the points z G Kt are lost from 
the domain of the map gt.) Theorem 1.2 (ii) states that, when Dc < D < Dc, two 
BES*^'^-' starting from different points < x < y < 00 can simultaneously return to 
the origin. In the complexificated version, all Z^{t) starting from z ^ Kt can arrive 
at the origin simultaneously {i.e., they are all swallowed). 

Osculation of the SLE path with M means that the SLE path has loops. Figure 
8 (a) shows the event that the SLE path makes a loop at time t > 0. The SLE 
hull Kt consists of the closed region encircled by the loop and the segment of the 
SLE path between the origin and the osculating point, and it is completely erased 
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9t 




Figure 7: (a) When 3/2 < L» < 2, the SLE^-^) path can osculate the real axis. The SLE hull 
is denoted by Kt- (b) The SLE hull Kt is swallowed. It means that all the points in Kt are 
simultaneously mapped to a single point —B{t) G M, which is the image of the tip of the SLE^^^ 
path 7(t). 

by the conformal transformation gt to EI as shown by Fig.8(b). Let < s < t and 
consider the map Qs, which is the solution of (1.20) at time s. Assume that 7(s) is 
located on the loop part of 7[0,t] as shown by Fig.8(a). The segment 7[0, s] of the 
SLE path is mapped by Qs to a part of M. Since 7(t) osculates a point in 7[0, s], its 
image gsi'jit)) should osculate the real axis M as shown by Fig.8(c). Since gj^ is 
uniquely determined from Qs, the above argument can be reversed. Then equivalence 
between osculation of the SLE path with M and self-intersection of the SLE path is 
concluded. In this intermediate phase < D < D^, 

SLE^^^ path 7 is self-intersecting, and 
|J;^ = H but 7[0, cx))nH7^e w.p.l. 

i>0 

[phase 3] When I < D < = 3/2, Theorem 1.2 (i) states for the Bessel flow that 
the ordering < is conserved for any < x < y. It implies that in this phase 
the SLE path should be a space-filling curve] 

7[0,oo) = e. 

(Otherwise, swallow of regions occurs, contradicting Theorem 1.2 (i).) 

Figure 9 summarizes the three phases of SLE paths. 

By complexification of Bessel flow, we can discuss flows on a two-dimensional plane H. 
By this procedure random curves (the SLE paths) are generated in the plane. The SLE 
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(b) 



(c) 



Figure 8: The event that the SLE^^^ path osculates M is equivalent with the event that the 
SLE^^) path makes a loop. 






(a) 



(b) 



(c) 



Figure 9: Schematic pictures of SLE^-'^) paths in (a) phase I {D > Dc = 2), (b) phase 2 
(Dc = 3/2 < L> < L»c = 2), and (c) phase 3 {l< D <Dc = 3/2). 



16 



paths are fractal curves and their Hausdorff dimensions ' are determined by Beffara 
[1] . We note that a reciprocity relation is found between D and ■* ; 

P-l)(dif)-l) = i D>D^ = '^-. (1.24) 

(In the phase 3, D < = 3/2, rfjf^ = 2.) 

Remark that the SLE map 7^ as well as the SLE path 7 are functionals of the BM. 
Therefore, we have statistical ensembles of random curves {7} in the probability space 
(O, J^, P). The important consequence from the facts that the BM is a strong Markov 
process with independent increments and gt gives a conformal transformation is that 
the statistics of {7} has a kind of stationary Markov property (called the domain Markov 
property) and conformal invariance with respect to transformation of the domain in which 
the SLE path 7 is defined [60]. 

The highlight of the theory of SLE would be that, if the value D is properly chosen, 
the statistics of {7} realizes that of the scaling limit of important statistical mechanics 
model exhibiting critical phenomena and fractal structures defined on an infinite discrete 
lattice. The following is a list of the correspondence between the SLE*^^^ paths with 
specified values of D and the scaling limits of models studied in statistical mechanics and 



fractal physics ^. 

gLE(3/2) ^ uniform spanning tree [62] 

SLE*^^/^) <^=^ critical percolation model 

(percolation exploration process) [80] 

SLE*^^-* <^==^ Gaussian free surface model (contour line) [78] 

SLE*^^/^) <^==^ critical Ising model (Ising interface) [81] 

SLE*^^/^) <^==^ self- avoiding walk [conjecture] 

SLE^^) «^=^ loop-erased random walk [62] 



^SLE^^) has a special property called the restriction property iff D = 5/2 (k = 8/3). It is well-know 
that the self- avoiding walk (SAW) model, which has been studied as a model for polymers, has this 
property. The conformal invariance of the scaling limit of SAW is, however, not yet proved. If it is 
proved, the equivalence in probability law between the scaling limit of SAW and the SLE^^/^) path will 
be concluded. 



17 



1.5 Dyson's BM model as multivariate extension of Bessel pro- 
cess 

Here we consider stochastic motion of two particles {Xi{t), X2{t)) in one dimension M 
satisfying the following SDEs, 

«.(t) = diS.W + f^^^^^^ t>0, (1.26) 

with the initial condition -^i(O) < -^2(0), where Bi{t) and B2{t) are independent one- 
dimensional standard BMs and /3 > is a 'couphng constant' of the two particles. 
The second terms in (1.25) represent the repulsive force acting between two particles, 
which is proportional to the inverse of distance X2{t) — Xi{t) of the two particles. 
Since it is a central force {i.e., depending only on distance, and thus symmetric for 
two particles), the 'center of mass' Xc{t) = {X2{t) + Xi[t))/2 is proportional to a BM; 

Xc(t) = B{t)/^/2 = B{t/2)^t > 0, where B{t) is a one-dimensional standard BM dif- 
ferent from Bi{t) and B2{t) and the symbol = denotes the equivalence in distribution. 
(Note that variance (quadratic variation) {dXc{t)f = {{dBi{t) + dB2{t))/2Y = dt/2, 
since {dBi{t)y — {dB2{t)Y — dt and dBi{t)dB2{t) — 0.) On the other hand, we can see 
that the relative coordinate defined by X,-{t) = (-^^2(0 — Xi{t))/V2 satisfies the SDE 

dX,{t) = dB{t) + ^^y t>0, 

where B{t) is a BM different from Bi{t), B2{t), B{t). It is nothing but the SDE for BES^^) 
with p + 

Dyson [16] introduced A'"-particle systems of interacting BMs in R as a solution X{t) — 
{Xi{t),X2{t), ...,XN{t)) of the following system of SDEs, 

dX,{t) = dB,{t) + ^ x(t)-X(ty ^ = 1'2,...,7V, (1.26) 

where {Bi{t)}f^ ^ are independent one-dimensional standard BMs and we define 

aij = inf{i > : Xi{t) ^ X^{t)), l<i <j <N, 

Tij = inf{t > aij : Xi{t) = Xj{t)}, l<i<j<N, 

T — min Ti-j. 
l<i<j<N ■' 

It is called Dyson's BM model with parameter /3 [65, 22]. As shown above, the N ~ 2 
case of Dyson's BM model is a coordinate transformation of the pair of a (time-change 
of) BM and BES*^^"*"^). In this sense, Dyson's BM model can be regarded as a multivariate 



18 



(multi-dimensional) extension of BES^^+^^Z^ > ^. In particular, we will characterize 
Dyson's BM model with ^ = 2 as an extension of the three-dimensional Bessel process, 
BES(=^), in Section 2. 



2 Two aspects of the Dyson model 

In this section, we study the special case of Dyson's BM model with parameter (3 = 2. 
We call this special case simply the Dyson model [48]. As shown above, the case (3 = 
2 corresponds to the case D = 3 oi Bessel process. In Sect. 1.2, we have shown that 
BES*^^^ has two aspects; (Aspect 1) as a radial coordinate of three-dimensional BM, and 
(Aspect 2) as a one-dimensional BM conditioned to stay positive. We show that the 
Dyson model inherits these two aspects from BES*^^^ [47]. 



2.1 The Dyson model as eigenvalue process 



Dyson introduced the process (1.26) with /3 = 1,2, and 4 as the eigenvalue processes of 
matrix- valued BMs in the Gaussian orthogonal ensemble (GOE), the Gaussian unitary 
ensemble (GUE), and the Gaussian symplectic ensemble (GSE) [16, 65, 22] 

For /3 = 2 with given e N, we prepare A"-tuples of one-dimensional standard 
BMs {-Bf/}^i, each of which starts from Xi e R, N{N - l)/2-tuples of pairs of BMs 



{Bij, Bij} 



l<i<j<N 



, all of which start from the origin, where the totally A^+2x A^(A^— 1)/2 



A^^ BMs are independent from each other. Then consider a,n N x N Hermitian-matrix- 
valued BM ^. 



( 



Bujt) - y/^Bi2{t) 
V2 



Bujt) + y/^Bujt) 



Bmjt) + ^^-^B^N{t) \ 

B2N{t) + ^^^B2N{t) 



BiN{t) - V^BiNit) B2N{t) - y/^B2N{t) 



\ 



^/2 



V2 



Km 



(2.1) 



^We can prove that r < oo for /3 < 1 and r = oo for /3 > 1 [73]. The critical value (3c = 1 corresponds 
to £)c = 2 of BES(^). 

^Precisely speaking, Dyson considered the Ornstein-Uhlenbeck processes of eigenvalues such that as 
stationary states they have the eigenvalue distributions of random matrices in GOE, GUE, and GSE. 
Here we consider matrix-valued BMs, so variances increase in proportion to time t>0. 

*In usual Gaussian random matrix ensembles, mean is assumed to be zero. The corresponding matrix- 
valued BM are then considered to be started from a zero matrix, i.e., Xi = 0,1 < i < N in (2.2). In random 
matrix theory, general case with non-zero means {i.e., Xi 0) is discussed with the terminology 'random 
matrices in an external source' (see, for example, [5]). From the view point of stochastic processes, 
imposing external sources to break symmetry of the system corresponds to changing initial state. 
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By this definition, the initial state of this EM is given by the diagonal matrix 



n^{Q)^dmg{xi,X2,...,XN). (2.2) 

We assume xi < X2 < ■ ■ ■ < x^. 

Remember that when we introduced BES*^^^ in Sect. 1.2, we considered the D-dimensional 
vector- valued BM, (1.4), by preparing D-tuples of independent one-dimensional standard 
BMs for its elements. Since the dimension of the space H(A^) of N xN Hermitian matrices 
is dimH(iV) = iV^, we need N"^ independent BMs for elements to describe a BM in this 
space H(iV). 

Corresponding to calculating an absolute value of B''^{t), by which BES^^^ was intro- 
duced as (1.5), here we calculate eigenvalues of T-i^if). For any i > 0, there is a family of 
N X N unitary matrices {U{t)} which diagonalize V-^{t), 

Uityn'^mit) = diag(Ai(i), . . . , Ajv(t)), t > 0. 

Let be the Weyl chamber of type A^v-i defined by 

W^ = {x^ (xi,X2,...,xn) eR^ : xi < X2 < ■■■ < xn}. 

If we impose the condition X{t) = (Ai(t), . . . , Ajv(t)) e W^,t > 0, U{t) is uniquely 
determined. _ 

For each t > 0, Xi{t),l < i < N are functional of {B^^{t), Bij{t), Bij}i<i^j<N, and 
then, again we can apply Ito's formula to take into account propagation of error correctly 
to derive SDEs for the eigenvalue process X{t),t > [7, 8, 41]. The result is the following, 

d\{t)=dB^^{t)+ Yl \.(t)-X.(ty ^^(O''^)' ^ = 1,2,..., AT, (2.3) 

where {B^"}f^^ arc independent BM difi^crcnt from the BMs used to define 'H^{t) by 
(2.1). It is indeed the /3 = 2 case of (1.26) as derived by Dyson (originally not by such a 
stochastic calculus but by applying the perturbation theory in quantum mechanics) [16]. 
Now the correspondence is summarized as follows. 

[Aspect 1] 

BES^^) <(=^ radial coordinate of 

D — vector- valued BM 

the Dyson model <^=^ eigenvalue process of 

with A^ particles N x N Hermitian-matrix-valued BM 

2.2 The Dyson model as noncolliding BM 

Here we try to extend the formula (1-11) for BES^^-* to multivariate versions. 
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First we consider a set of two operations, identity (cJi = id) and reflection ((72 = ref), 
such that for x e R, ai{x) — x and (72 (x) = —x, and signatures are given as sgn((7i) = 1 
and sgn((72) — —1, respectively. Then we have 

pf'{y\x)= Yl sgn{a)pt{y\a{x)), t>Q, x,y,eR. (2.4) 

iTG{id,ref} 

Next we consider a set of all permutations of N indices {1,2,..., N}, which is denoted 
by <Sjv, and put the following multivariate function 



N 

V sgn{a)Y\pt{ya(i)\xi) = ^ det Pt{yi\xj) (2.5) 



of cc = {xi, . . . ,xn) e W;^ and x = [yi, . . . , y^) G with a parameter t > 0. Following 
the argument by Karlin and McGregor [34], we can prove that this determinant gives 
the transition probability density with duration t from the state x to the state y of A^- 
dimensional absorbing BM, B^{t) = (fif (t), . . . , fi^^(t)), in a domain W;^, in which 
absorbing walls are put at the boundaries of W^. Since the boundaries of are the 
hyperplanes Xi = Xj,l<i<j<N, the Brownian particle B^{t) is annihilated, when any 
coincidence of the values of coordinates of (t) occurs. The 'survival probability' that 
the BM is not yet absorbed at the boundary is a monotonically decreasing function of time. 
If we regard the i-th coordinate B^' (t) as the position of i-th particle on M, 1 < i < A^, the 
state X e is considered to represent a configuration of N particles on M such that a 
strict ordering Xi < X2 < ■ ■ ■ < xn of positions is maintained, while the state absorbed at 
the boundary, x e dW^, is a configuration in which collision occurs between some pair of 
neighboring pairs of particles; 1 <^ i < iV — 1, s.t. Xi — Xi+i. Since if such collision occurs, 
the process is totally annihilated, this many-particle system is called vicious walker model 
[20, 57, 32, 38, 33, 12]. 

As already noted below Eq.(l.ll), hi{x) = a; is a harmonic function in (0, 00) condi- 
tioned 0(0) = 0. Similarly, if we consider a harmonic function of N variables 

i=l * 

conditioned hN{x) — 0,^x & dW^, we will have the following product of differences 

hN{x)= Yl i^j-^i)^ (2-7) 

i<i<j<N 

which is identified with the Vandermonde determinant deti^ij^N[xi~^]. 
Combining above consideration, we put the following function 

pf(y|a;) = ^^ det L(|/,|a;,)l, t>0, x,y eW^. (2.8) 
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We can show that the factor hN{y)/hN{x) provides an exact renormahzation of the vicious 
walker model (the absorbing BM in W^) to compensate any decay of total mass of the 
process by collision (by absorption at 9W;^), and that (2.8) gives the transition probability 
density function for the A^-particle system of one-dimensional BMs conditioned never to 
collide with each other forever (which we simply call the noncolliding BM) [25, 38, 39]. 

Moreover, by using the harmonicity (2.6), we can confirm that (2.8) satisfies the 
following partial differential equation, 

l<i^j<N * ^ * 

with the initial condition Pq {y\x) = 5{y — x) = YliLi ^iVi ~ ^i) [38]- It can be regarded 
as the backward Kolmogorov equation of the stochastic process with N particles, X{t) = 
{Xi{t), . . . , Xiv(i)), which solves the system of SDEs, 

dX,{t) = dBt^{t)+ X(t)-X(ty ^^^0'^)' ^ = 1>2,...,7V. (2.10) 

Eq.(2.10) is identified with the /3 — 2 case of (1.26). Then the equivalence between 
the Dyson model and the noncoUiding BM is proved. 
The result is summarized as follows. 
[Aspect 2] 

BES^^-* <^=^ /i-transform of absorbing BM in (0, oo) 

-^=^ BM conditioned to stay positive 

the Dyson model <(=^ /i-transform of absorbing BM in 

noncolliding BM 

The fact that BES^^^ and the Dyson model have two aspects imphes useful relation 
between projection from higher dimensional spaces and restriction by imposing condi- 
tions. For matrix-valued processes, projection is performed by integration over irrelevant 
components, and noncolliding conditions are generally expressed by the Karlin-McGregor 
determinants. The processes discussed here are temporally homogeneous ones, but we can 
also discuss two aspects of temporally inhomogeneous processes. Actually, we have shown 
that, from the fact that the temporally inhomogeneous version of noncolliding BM has 
the two aspects, the Harish-Chandra-Itzykson-Zuber integral formula [26, 30] is derived. 



/ 



r 1 ^ C a^^ r 1 

dU exp <^ -— Tr(Aa; - U*KyUf \ = \ det pt{y,\x,) , (2.11) 

[ 2cr^ ' J hN{x)hN{y) i<i,J<N L ' "1 



'U(iV) 

(7^ > 0, where dU denotes the Haar measure of U(A'^) normalized as /u(;v) dU = 1, Ax — 

diag(a;i, . . . , x^v), Ay = diag(yi, . . . .Vn) with x,y e W%, and Cn = (27r)^/^ Hil^ r(i) 
[40] 9. 



^We can apply the present argument also for processes associated with Weyl chambers of other types. 
See [49, 41]. 
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3 Determinantal processes and entire functions 



3.1 Aspect 2 of the Dyson model 

As Aspect 2, the Dyson model is constructed as the /i-transform of the absorbing BM in 
W^. Therefore, at any positive time t > the configuration is given as an element of 

x{t) = (Xi(t), X2(t), . . . , e w^, t> 0, (3.1) 

and there is no multiple point at which coincidence of particle positions, Xi{t) = Xj[t),i ^ 
j, occurs. That is, the Dyson model is equivalent with the noncolliding BM. We can 
consider the Dyson model, however, starting from initial configurations with multiple 
points. In order to describe configurations with multiple points, we represent each particle 
configuration by a sum of delta measures in the form 

C{-) = Y.^M t>0 (3.2) 

with a sequence of points in R, a; = (a^i)iGn, where I is a countable index set. Here for 
y e R, Sy{-) denotes the delta measure such that Sy{x) = 1 for x = y and Sy{x) = 
otherwise. Then, for (3.2) and A C M, ^{A) = f^^{dx) = J^iei-xteA'^ ~ • ^ 
If the total number of particles N is finite, I = {1,2,..., N}, but we would like to also 
consider the cases with N — oo. We call measures of the form (3.2) satisfying the condition 
< oo for any compact subset K G M. nonnegative integer- valued Radon measures 
on R and write the space of them as DJl. The set of configurations without multiple point 
is denoted by TIq = G 9Jt : ^{{x}) < l,^x G R}. There is a trivial correspondence 
between and 2Jlo- 

First we assume ^ = I^iei ^ 37to, C(I^) = N eN and consider the Dyson model as 
an QJto-valued diffusion process, 

E{t,-) = J2^xut){-), t>Q, (3.3) 

starting from the initial configuration ^, where X{t) = {Xi{t), • • • , Xjv(i)) is the solution 
of (2.10) under the initial configuration x — {xi, . . . , xn) G W;^. We write the process as 
(S(t),P«) and express the expectation with respect to the probability law of the Dyson 
model by IE^[-]. We introduce a filtration {J^(t)}fe[o,oo) on the space of continuous paths 
C([0, oo) — )■ DJl) defined by J-'(t) = a{E{s), s G [0, t]), where a denotes the sigma field. 

Then we introduce a sequence of independent one-dimensional standard BMs, B''^{t) = 
{Bf'(t))i^i, t >0 and write the expectation with respect to them as E^[ • ]. 

Let liyyA(a;) = 1 if a; G and 1^a{x) = otherwise. Then Aspect 2 of the Dyson 
model is expressed by the following equality; for any < t < T < oo, any symmetric 
function g on R^, 

E^[g{Xm = E 



\hN{Bm 
hN{x) 



(3.4) 



23 



where we have assumed the relations ^ = X^jgj G 2?to; = N e N,x = {xi, . . . ,xn) G 
and (3.3). Note that the indicator I^a in the RHS annihilates the BM, B{t), 

if it hits any boundary of the Weyl chamber, dW^, and the factor \h^r (B {t)) \ / (x) per- 
forms the /i-transform of measure for Brownian paths. That is, the RHS of (3.4) indeed 
gives the expectation of g with respect to the process obtained by the /i-transform of the 
absorbing BM in W^. 

If we apply the Karlin-McGregor determinantal formula (see (2.5)), 

(RHS) = E' 

where we recall the definition of determinant and let a{B{t)) = (-Bcr(i), . . . ,Bcr(N)),o' G 
Sn- Since hN{x) is a product of differences, ^^^^^ sgn((j)liyyA ((T(S(i)))|/iAr(-B(i))| — 
YIic^Sn '^wj^{'^{^{'t)))hN{B{t)) — hN{B{t)), and then (3.4) is simply written as 

In the LHS of (3.5) we should note that the Dyson model is an interacting particle 
system such that between any pair of particles a long-range repulsive force acts. Strength 
of the two-body repulsive force is exactly proportional to the inverse of distance between 
two particles and thus it diverges as the distance goes to zero. By this strong repulsion, 
any coUision of particles is prevented. On the other hand, in the RHS of (3.5), independent 
BMs are considered. When we calculate the expectation of a function g of them, however, 
we have to put extra weight hM{B{t)) / h]\j{x) to their paths. Since if \Bj{t) — Bi{t)\ 
for any i ^ j, this weight becomes zero, again any collision of particle is prevented. An 
important point of the Karlin-McGregor formula is that this weight for paths is signed, 
i.e., it can be positive and negative. Therefore, all particle configurations realized by 
intersections of paths in the 1-1-1 spatio-temporal plane are completely cancelled. 

3.2 Complex BM representation 

Now we consider complexification of the expression (3.5) [48]. For each B^^{t),i G I, we 
introduce an independent one-dimensional BM starting from the origin, Bi{t), and define 
a complex BM as Zi{t) = Bi(t) + \^ —lBi(t) , i G I. If we write the expectation with respect 
to {Bi{t)}i(,i as E[ • ] and define = E, we can confirm that the RHS of (3.5) can 
be replaced by 

(3.6) 

where Z{t) = {Zi{t))iei. 



J2 sgn(a)lwA(a(S(t)))^^(B(t))M^ 



9iB{t)) 



hN{B{t)) 
/iiv(cc) 



(3.5) 



E 



X 



gm)) 



hN{Z{t)) 
hN{x) 
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A key lemma of our theory [48] is the following identity; for a; = (xi, . . . , xat) e W^, 



hNix) 



det 

l<i,j<N 



where 



= n 

Xgsupp i;n{u}'^ 



z — u 
X — u 



(3.7) 



for ^ = e SOTo and supp ^ = G M : i{{x}) > 0}. The function ^^(z) has 

an expression of the Weierstrass canonical product with genus zero. Then, it is an entire 
function with zeros at supp ^ Pi {uY (see, for example, [64, 69]). 

If we apply this identity to (3.6), we have quantities ^^*{Zj(t)),i, j G I, which are 
conformal transforms of independent complex BMs, Zj{t),j G I. Since complex BM is 
conformal invariant, each $|*(Zj(t)) is a time change of a complex BM, Zj{-). Then the 
average is conserved. 



E^[$^(Zj(i))] = E^[$^(Z^-(T))], <^ t<T <oo, 



(3.8) 



that is, {$^(Zj(t))}jg]i are independent conformal local martingales (see, for example. 
Section V.2 of [72] ). 

Let < t < T < oo. Then for any J-'(t)-measurable function F in the continuous path 
space C([0, oo) Wl), we have the equality 



E«[F(H(-))] = E 



X 




$-(Z,(T)) 



(3.9) 



Now it is claimed that any observables of the Dyson model is calculated by a system of 
independent complex BMs, whose paths are weighted by a multivariate complex function 
detij£i[$|*(Zj(T'))], which is a conformal local martingale. We call (3.9) the 'complex 
BM representation' of the Dyson model [48]. 



3.3 Determinant al process with an infinite number of particles 

For a configuration ^ G 9Jl, we write the restriction of configuration in A C R as (^ fl 

ie^^■.x^eA^xii^)^ ^ shift of configuration by m G M as t„^(-) — ^jgj (^a;j+u(')' a 
square of configuration as ^^^^(O = X]iei^^^(')' respectively. Let Co(lR) be the set of all 
continuous real-valued functions with compact supports. 

For any integer M G N, a sequence of times t = {ti, ^2, • • • , ^m} with < ti < • ■ ■ < 
tM < T < oo, and a sequence of functions / = (/j^, ft^, . . . , ftj^) G Co(R)^, the moment 
generating function of multitime distribution of (S(t),P^) is defined by 



CXp<^ / ftm{x)-itm,dx) 



(3.10) 
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Prom the fact that {^'^{Zj{t))}j(zi are independent conformal local martingale with 
the property (3.8), 



E 



X 



$|'(z,(r)) 



= E 



X 



$-(Z,(0)) 



Moreover, we can see that the following property holds for the determinantal weight in 
the complex BM representation (3.9). Let f be a subset of the index set I and assume 
that a function / depends on Bi{t),i e F, but does not on Bj{t),j e < i < T < oo. 
Then 



E 



X 



/({5,(t)},eF)det[$-(Z,(T))] 



E 



X' 



i,jei' ^ 



(3.11) 



where E*'^'[-] denotes the expectation with respect to {i?^''(t)}jgF. Let ^ e TXq and 



Jr Jr 

-l{s > t)ps-t{x\y), (s, t) e (0, oof, {x, y) e M^ 



(3.12) 



where l(a;) is the indicator function of a condition ui] l(a;) = 1 if a; is satisfied and 
1(0;) = otherwise. Then, using this 'reducibility' (3.11), we have proved that (3.10) is 
given by a Fredholm determinant 



Det 

(s,^)et^ 
(a;,2/)eR2 



^stK^ -y) + IK^("S, X] t, y)xt{y) 



(3.13) 



where Xtm(') = e^'^"^^'^ — 1, 1 < m < M. We call a correlation kernel. 
For L > 0, a > and ^ e 9Jl we put 



Odx) 

L,L]\{0} X 



C{dx] 

L,L]\{0} 1^1 



1/a 



and M(0 = lim^^oo M(^, L), M„(0 = lim^^oo M«(^, L), if the limits finitely exist. We 
have introduced the following conditions for initial configurations ^ E DJt [45]: 

(C.l) there exists Cq > such that \M{^, L)\ < Co, L > 0, 

(C.2) (i) there exist a G (1, 2) and Ci > such that M^iO < Ci, 
(ii) there exist /3 > and C2 > such that 



Mi(r_„2^<2>) < C2(max{|a|, 1})"'^ Va e supp 
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It was shown that, if ^ G TXq satisfies the conditions (C.l) and (C.2), then for a e R and 
zeC, = Jhn $|n[a-L,a+L](^) finitely exists, and 

< Cexp|c(|a|^ + l^nj i ^ a e supp ^, z E C, 

for some c, C > and 9 G (max{a, (2 — /3)},2), which are determined by the constants 
Co,Ci,C2 and the indices a,/3 in the conditions [45]. Then even if ^(M) = oo, under 
the conditions (C.l) and (C.2), given by (3.12) is well-defined as a correlation kernel 
and dynamics of the Dyson model with an infinite number of particles (S(t),P^) exists 
[45]. We note that in the case that ^ G OJIq satisfies the conditions (C.l) and (C.2) with 
constants Co,Ci,C2 and indices a and (3, then ^ fl [— L,L],VL > does as well. Then 
we can obtain the convergence of moment generating functions "^^^ ^'^'[/] ^|[/] ^ 
L — >■ oo, which implies the convergence of the probability measures P^'^[~^'^l — >■ in 
L — )■ oo in the sense of finite dimensional distributions [45] . 

By definition of Fredholm determinant, the moment generating function (3.13) can be 
expanded with respect to Xt^(-), 1 <m < M, as 



Nru>0, 

l<m<M 



M 



Nrr 



n 

L In 



"Nr, 



(m) 
Nm 



m=l 



1=1 



X. 



(m) 



{M)\ 
Nm )■ 



with 



(1). .4. 

"^Ni ) ■ ■ ■ ) '^Mi a? 



{My 

Nm , 



det 

l<i<Nm,l<j<N„, 

l<m,n<M 



where x"^^ denotes {x^i^\ . . . , I) and dx'^^l = n-L"[ dx^'"' , 1 < m < M. The functions 



p^'s are called multitime correlation functions, and ^|[/] can be regarded as a generating 
function of them. In general, when the moment generating function for the multitime dis- 
tribution is given by a Fredholm determinant, all the spatio-temporal correlation functions 
are given by determinants of matrices, whose entries are special values of a continuous 
function K^, and then the process is said to be determinantal [68, 37, 43]. (Therefore, 
is called a correlation kernel.) The results by Eynard and Mehta reported in [17] 
for a multi-layer matrix model can be regarded as the theorem that the Dyson model is 
determinantal for the special initial configuration ^ = N6q, i.e., all particles are put at 
the origin, for any N e N. The correlation kernel is expressed in this case by using the 
Hermite orthogonal polynomials [67]. The present author and H. Tanemura proved that, 
for any fixed initial configuration ,^ G OJt with .^(M) G N, the Dyson model (5(i),P^) is 
determinantal, in which the correlation kernel is given by 

K^{s,x;t,y) = / dzps{x\z) [ dwpt{-\f-ly\w) 

(l - ^lEh^] _ i(, > t)ps-t(x\y), (3.14) 
\ X — z J 



n 



x'Gsupp ^ 
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where r(^) is a closed contour on the complex plane C encircling the points in supp ^ on 
the real line R once in the positive direction [45]. When ^ e QJIq, (3.14) becomes (3.12) 
by performing the Cauchy integrals. 

We note that = X]te£ ® S(^) is a determinantal point process (or Fermion point 
process) on the spatio-temporal field t x M with an operator /C given by ICf{s,x) = 
^jg^ /j^(iy]K(s,a;;t,|/)/(t,|/) for /(t, ■) G Co(M),t G t. When /C is symmetric, Soshnikov 
[82] and Shirai and Takahashi [79] gave sufficient conditions for to be a correlation 
kernel of a determinantal point process (see also [27]). Such conditions are not known 
for asymmetric cases. The present correlation kernels (3.12) and (3.14) are asymmetric 
cases, K.^{s,x;t,y) ^ ]K^(t, s, a;) by the second terms — l(s > t)ps„((x|?/). Such form of 
asymmetric correlation kernels is said to be of the Eynard-Mehta type [4, 48]. 

Prom the view point of statistical physics, such asymmetry is useful to describe 
nonequilibrium systems developing in time. In order to demonstrate it, we have stud- 
ied the following relaxation phenomenon of the Dyson model with an infinite number of 
particles. 

We consider the configuration in which every point of the integers Z is occupied by 
one particle. 

See Fig.lO. It can be confirmed that satisfies the conditions (C.l) and (C.2) and thus 
the Dyson model starting from (S(t),P^ ), is well-defined as a determinantal process 
with an infinite number of particles. 

As a matter of fact, we have shown that the correlation kernel is given by [45] 

K^^{s, x; t, y) = Ksin(s, x; t, y) 

+— / dke'''^'-'^l''+'^^^y-''^U^{x-^/-^ks,2^:^/-^s)-l\ 
27r JifcK^ I- J 



{s,x]t,y) 



neZ\{0} 



27r\ 



TTu{{y — x) — 27rv — Isw} 



{s,t) e [0,oo)^, {x,y) e M^, where 



Ksin{s,x;t,y) 



2n 



- f dk efc'(*--)/2+^My--) _ > t)p, 

[ due''^''^^^-'^/^cos{7iu{y-x)} ii t 
Jo 

Ksin{x,y) 

POO 

- du e'^'"'(*-^)/2 cos{TTu{y - x)} iit<s, 
Ji 



x\y) 



> s 



iit^s 
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Figure 10: Consider tlie Dyson model starting from the configuration in which every point of 
the integers Z is occupied by one particle. This nonequilibrium determinantal process shows a 
relaxation phenomenon to the stationary state fisin- 

with 

KUx,y) = ^ [ dke^'^^y-^^ = '^"^i^y-^)} ^ e M, (3.15) 

and "i^a is a version of the Jacobi theta function defined by 

^?3(w,r) = ^e2-^^"+"^™', Qt>0. (3.16) 

By this exphcit expression, we can see that 

hm K^^{u + s, x;u + t,y) = Ksin(s, x; t, y). (3-17) 

U—>OD 

The correlation kernel Ksin(s, x; t, y) surviving in the long-term limit is called the extended 
sine kernel It is symmetric, Ksin(s, x; t, y) = Ksin(t, y; s, x), and the determinantal process 
with the correlation kernel Ksin(s, x; t, y) is an equilibrium dynamics. It is time-reversal 
with respect to the determinantal point process /isim in which any spatial correlation 
function is given by a determinant with the correlation kernel (3.15) called the sine kernel. 
This stationary state iisin is a scaling limit (called the bulk-scaling-limit) of the eigenvalue 
distribution of random matrices in GUE [65, 22]. See [83, 70, 67, 43]. 

The theory of entire functions discusses the relations between the growth of an entire 
function and the distribution of its zeros [64, 69]. Here we set distributions of zeros 
as initial configurations satisfying some conditions and control the behavior of particles 
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at infinity to realize nonequilibrium dynamics of long-rang interacting infinite-particle 
systems. Systematic study of determinantal processes with infinite number of particles 
exhibiting relaxation phenomena is now in progress [44, 45, 46, 48] 

4 Related Topics 

At the end of this manuscript, I briefly introduce related topics, which we are interested 
in. 

4.1 Extreme value distributions of noncoUiding diffusion pro- 
cesses 

As explained in Section 1.2, BES*^^^ is the conditional BM to stay positive. When we 
impose an additional condition such that it starts from the origin at time t — and 
return to the origin at time t — 1, the process is called the three-dimensional Bessel 
bridge with duration 1, which is here denoted hy Y{t),0 < t < 1. 
Here we consider the maximum of Y{t), 

Hi = max Y(t). 

o<t<i ^ ' 

See Fig. 11 We can show that the distribution of H\ is described as [52] 

P(i^i <h) = J]i/2(\/2n/i)e-2'^'"', (4.1) 



where Hi[x) is the i-th Hermite polynomial 

-l)^(2a 



H^i^)=''■T. \v■-tM ^ {0,1,2,...} (4.2) 



j=0 



with [a]=the greatest integer that is not greater than a G M. By the equation P(i?i < 
h) = Jq pi{u)du, the probability density function pi{u) of Hi is defined, and s-th moment 
of the random variable Hi is calculated by 



POO 

E[Hl] = / dhh'piis). 
Jo 



^°It will be interesting to discuss intrinsic relations between the above mentioned relaxation phe- 
nomenon to /Usin and the interpolation/sampling theorem of Whittaker and others [9], which represents a 
function / in the cardinal sampling series 



TT z — n 
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Time 



Figure 11: Let Hi be the maximum of the three-dimensional Bessel bridge with duration 1. 
The expectation E[iJf ] is expressed by using the value of Riemann zeta function at s G C. 



As discussed by Biane, Pitman, and Yor [2], if we set E[iff] = 2(7r/2)*/^^(s), the following 
equality is established, 

as) = ls{s-l)n-^/'T{s/2)as) 

= ^ + ^s(s-l)J^ du(u'/'-' + u^'-'y^-'){UO,V^u)-l), (4.3) 

where r{z) is the gamma function, ({z) is the Riemann zeta function, 



oo _ 



(4.4) 



n=l 



and ■i?3(v,r) is given by (3.16). See also Chapter 11 in [87]. 

Then we consider an A^-tuples of three-dimensional Bessel bridges conditioned never 
to collide with each other, Y^^\t) = {Yl^\t), Y}f\t)), 

o<F(^)<rf)<---<rf\ 0<t<l, 

with Y^^\0) = Y^^\l) = 0. It is found that (4.1) is generahzed for 

— max yi^^ 
o<t<i ■'^ 



as [76, 52] 



P{Hn < h) 



:-i 



det 



2^' n^=ir(2i) i<ij<^ 



2{j+j-l) 



{\/2nh)t 



.neZ 



(4.5) 



Extensive study of extreme value distributions of noncolliding processes has been reported 
in [24, 36, 18, 76, 19, 3, 66, 71, 31]. See also [53, 86, 42, 58]. 



31 



Recently, Forrester, Majumdar, and Schehr clarified an equality between the distri- 
bution function P(ifjv < h) and the partition function of the 2d Yang-Mills theory on a 
sphere with a gauge group Sp(2A^) [23]. Moreover, by using this equivalence, they proved 
that in a scaling limit of P(-ffjv < h), the Tracy- Widom distribution of GOE of random 
matrices [84, 85] is derived. 



4.2 Characteristic polynomials of random matrices 

Let "H be an X Hermitian random matrix in the GUE. The characteristic polynomial 
of a variable a is then given by 

P{a) ^det{al -n), 

where / is the N x N unit matrix. In the connection with the Riemann zeta function 
(4.4), statistical property of P{a) has been studied [50, 51, 28, 6]. Here we consider the 
ensemble average of m-product of characteristic polynomials 

MGUE(m,a;7V,a2) = /j]p(Q;„)\ ^ {T[\[{ocn - K)) , (4.6) 

\n=l / GUE(Ar,(T2) \n=l i=l / GVE{N,(j'^) 

where {■)gue{n<j^) denotes the ensemble average in the GUE with variance a'^ oi N x 
N Hermitian matrices {K}, whose eigenvalues are {{K)iLi}- The probability density 
function of eigenvalues of T-L in GUE with variance is given by 

I^nAO = exp (--^j hNixf, (4.7) 

e = EliSx. em,x,<X2<---<XN, where Cn = (27r)^/2 Jlti r(i), \x\' = EZi^I 
and /i7v(a?) is given by (2.7). Then if the expectation of a measurable function F of a 
random variable S e M with respect to (4.7) is written as 

with setting ^ = ^^j^ (5a;., a; = {xi, . . . ,xn), where dx — YliLi dxi, (4.6) is written as 



MGVEim,a;N,a^) = E^,^2 



,n=l XeS 



In Section 3.3, we showed that the Dyson model (S(t),P^) is a determinantal process 
with the correlation kernel (3.12) for any fixed initial configuration ^ G 9Jto if ^{^) = 
N eN. Here we consider the situation such that the initial configuration ^ is distributed 
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according to (4.7). Note that by the term hNixf in (4.7), the GUE eigenvalue distribution 
is in OJto w.p.l. 

By using Aspect 2 of the Dyson model, we can prove that this process, denoted by 
(S(t), P^'^''^^), is equivalent with the time shift t — )> t + o"^ of the Dyson model starting from 
the configuration A^^o {i.e., all particles are put at the origin). That is, the equality 



(S(i),P'^^.-^) = (S(i + (72),P^^°) 



(4.8) 



holds for arbitrary o"^ > in the sense of finite dimensional distribution [35]. 

This equivalence is highly nontrivial, since even from its special consequence, the 
following determinantal expression is derived for the ensemble average of product of char- 
acteristic polynomials; For any A^, n e N, a = (ai, 0:2, ■ ■ ■ , Q!2n) G > 0, 



MGUE(2n,a;Ar, a^) 
1 



n{2N+n) 



X det 

l<ij<n 



Gtq Gt 



n+j 



hn{ai, ■ ■ ■ , a„)/i„(a„+i, ■ ■ ■ , a2„) 

0-2) HN+nian+j/V^) 



H 



N+n 



(4.9) 



where 



i=2 



{N + n-i)\ 
{N + n-l)V 



and Hi{x) is the i-th Hermite polynomial given by (4.2). 

Moreover, in order to simplify the above expression, we can use the following iden- 
tity, which was recently given by Ishikawa et al. [29] as a generalization of the Cauchy 
determinant 



det 



1 



hn{x)hn{y) 



l<id<n \Xi + yj J ^^=l n?=l i^i + yj)' 

For n > 2, a; = {xi, . . . , x„), y = (yi, . . . , y„), a = (oi, . . . , a„), 6 = (61, . . . , 6„) e C", 
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(4.10) 



Then we have the expression 

MGUE(2n, a; N, a^) 



1 



h2n{Oi) l<j,i<2n 



HN+i-i{o:j; 



(4.11) 
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where 

This determinantal expression (4.11) is also obtained from the general formula given by 
Brezin and Hikami as Eq.(14) in [6]. See [13] for recent development of this topic. 

4.3 Fomin's determinant for loop-erased random walks and its 
scaling limit 

We consider a network F = (V, E, W), where V — {vi} and E = {cj} arc sets of vertices 
and of edges of an undirected planar lattice, respectively, and W = {w{e)}eeE is a set of 
the weight functions of edges. For a, 6 e let tt be a walk given by 

ei 62 63 e™ , 
TT : a — Vq Vi V2 ■ ■ ■ ^ Vm — 

where the length of walk is |7r| = m e N and, for each < i < m — 1, Vi and Vi+i are 
nearest-neighboring vertices in V and ei E E is the edge connecting these two vertices. 
The weight of tt is given by w{7r) — YllLi'^i^i)- For any two vertices of a, 6 e V, the 
Green's function of walks {tt : a — > 6} is defined by 

w{a,b) = j2 E ^w- 

m 7r:a— >fe,|7r|=m 

The matrix W = {W{a, b))a^bev is called the walk matrix of the network F. 

The loop-erased part of tt, denoted by LE(7r), is defined recursively as follows. If 
TT does not have self-intersections, that is, all vertices Vi,0 < i < m are distinct, then 
LE(7r) = TT. Otherwise, set LE(7r) = LE(7r'), where n' is obtained by removing the first 
loop it makes. The loop-erasing operator LE maps arbitrary walks to self- avoiding walks 
(SAWs). Note that the map is many-to-one. For each SAW, the weight w{() is given 

by 

w{c)= Yl ^w- (4.12) 

7r:LE(7r)=C 

We consider the statistical ensemble of SAWs with the weight (4.12) and call it loop-erased 
random walks (LERWs) [61]. 

Assume that A = {oi, a2, . . . , aj^} C V and B — {bi, b2, ■ ■ ■ , bj^} C V are chosen so 
that any walk from to bj intersects any walk from aj/, i' > i, to bj',j' < j. The weight of 
A^-tuples of independent walks ai — ^ 6i, . . . , ajv is given by the product of weights 

YliLi ""^("""i)- Then we consider iV-tuples of walks (tti, 112, ... , tttv) conditioned so that, for 
any 1 < i < j < N, the walk ttj has no common vertices with the loop-erased part of tt^; 

LE(7ri) n TTj- = 0, l<i<j<N. (4.13) 
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Figure 12: The situation LE(7rj) n vrs = 0, j = 1, 2 is illustrated in a planar domain D, where 
A = (01,02,03) and B = (61,62,^3) are ah boundary points of dD. In this figure, LE(7ri) 
and LE(7r2) denoted by solid curves are the loop-erased parts of the walks 711:01—7- hi and 
7r2 : 02 — 62) respectively. The third walk vrs : 03 — )• 63 can be self-intersecting, but it does not 
intersect with LE(7ri) nor LE(7r2). 

See Fig. 12. By definition, LE(7rj) is a part of vTj, and thus nonintersection of any pair of 
loop-erased parts is concluded from (4.13); 

LE(7r,) n LE(7rj) = 0, l<i<] <N. 

Fomin proved that total weight of A^-tuples of walks satisfying such a version of nonin- 
tersection condition is given by the minor of walk matrix, detiWA^s) = deta,^A,bi^BiW {a, b)) 
[21]. This minor is called Fomin's determinant and Fomin's formula is expressed by the 
equality [21, 61] 

N 

det{WA,B)= Yl ll^i^k). (4.14) 

LE(7ri)n7rj=0,i<j fc=l 

Kozdron and Lawler [55] consider continuum limit (the diffusion scaling limit) of 
Fomin's determinantal system of loop-erased random walks in the complex plane C, 
where the initial and the final points A = {aj} and B = {6,} of paths can be put on 
the boundaries of the domains dD. By the diffusion scaling limit each random walk will 
converge to a path of complex BM. We should note that, however, the characteristics of 
BM look more similar to those of a surface than those of a curve. It implies that the 
Brownian path has loops on every scale and then the loop-erasing procedure mentioned 
above does not make sense for BM in the plane, since we can not decide which loop is 
the first one. Kozdron and Lawler proved explicitly, however, that the continuum limit 
of Fomin's determinant of the Green's functions of random walks converges to that of the 
Green's functions of Brownian motions [55]. This will enable us to discuss nonintersecting 
systems of loop-erased Brownian paths in the sense of Fomin (4.13). Moreover, Kozdron 
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[54] showed that 2x2 Fomin's determinant representing the event LE(/9i) n/52 — for two 
complex Brownian paths (^1,^2) is proportional to the probability that 7sle(3) H ^ = 0, 
where 7sle(3) and (3 denote the SLE^^^ path and a complex Brownian path, respectively. 
On the other hand, Lawler and Werner gave a correct way to add 'Brownian loops' to an 
SLE*^^) path to obtain a complex Brownian path [63]. These results imply that the scaling 
limit of loop-erased part of complex Brownian path is described by the SLE^^^ path, as 
announced in Section 1.4. 

Setting a sequence of chambers in a planar domain, Sato and the present author ob- 
serve the first passage points at which A^-tuples of complex Brownian paths . . . , /J^r) 
first enter each chamber, under the condition that the loop-erased parts (LE(/3i), . . . , LE(/3jv)) 
make a nonintersecting system in the domain in the sense of Fomin (4.13) [75]. It is proved 
that the system of first passage points is a determinantal point process in the planar do- 
main, in which the correlation kernel is of Eynard-Mehta type [75] . Interpretation of this 
result in terms of 'mutually avoiding SLE paths' [56, 15] will be an interesting future 
problem. 
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